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Finite-top-mass effects in NNLO Higgs production 
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We construct an accurate approximation to the exact NNLO cross section for Higgs production in gluon-gluon 
fusion by matching the dominant finite top mass corrections recently computed by us to the known result in 
the infinite mass limit. The ensuing corrections to the partonic cross section are very large when the center 
of mass energy of the partonic collision is much larger than the Higgs mass, but lead to a moderate correction 
at the percent level to the total Higgs production cross section at the LHC. Our computation thus reduces the 
uncertainty related to these corrections at the LHC from the percent to the per mille level. 



The search for the Higgs boson is one of major 
tasks of the forthcoming experiments at the Large 
Hadron Collider (LHC) at CERN. The theoretical 
and experimental effort which has been put into 
Higgs studies for LHC phenomenology is remark- 
able. In particular, the determination of higher- 
order corrections in perturbative QCD has been 
widely investigated. The dominant Higgs produc- 
tion mechanism in the Standard Model is gluon- 
gluon fusion through a top loop. The hadronic 
cross section can be obtained by convolution of 
the partonic cross section with parton distribu- 
tions Uix, 
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^ dxi dx 



Xi 



X2 



Cij as; 



Th 

X1X2 ' 



, yt fi (a;i, m|^) /, {x2,mjj) (1) 



where (7o{yt) is the partonic Born cross sec- 
tion [T] and the dimensionless variables r/j and yt 
parametrise the hadronic center-of-mass energy 
and the dependence on the top mass: 



Th = 



yt 



m 



(2) 



H 



The dimensionless coefficient function 



C{as;T,yt) contains the QCD corrections. The 
NLO contribution to it was computed in [5] and 
recently confirmed in [3]. The dominant NLO 
correction comes from the radiation of soft glu- 
ons, which cannot resolve the quark loop in the 
ggH coupling. Therefore, at least at the inclu- 
sive level, the approximation to the exact NLO 
result obtained [4I5| by taking the limit mt — > 00 
turns out to be very accurate. This approxi- 
mation considerably simplifies the calculation 
because the ggH coupling becomes pointlike and 
the corresponding Fcynman diagrams have one 
less loop. 

Recently, the NNLO contribution to 
C{as]T,yt) has been computed in the — » cxd 
limit [5]. The NNLO result appears to be per- 
turbatively quite stable and it should provide 
a good approximation to the yet unknown ex- 
act result; it has been widely used for precision 
phenomenology at the LHC l^. However, the in- 
finite mt approximation fails in the limit of large 
partonic center-of-mass energy or equivalently 
T ^ 0. This is due to the fact that the high 
energy behaviour of the partonic cross section 
is completely different according to whether the 
ggH coupling is pointlike or goes through a quark 
loop, because in the latter case the quark loop 
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Figure 1. The NLO coefficient function for 
rriH — 130 GeV. The red curve corresponds to 
the exact case, the black one to mt — > oo and 
the blue one to the approximation Eq. ([6]) with 
To = 0.057 and uj = 1/100. 

effectively provides a form factor which softens 
the interaction. Indeed, the coefficient function 
behaves respectively as 



°° 1„2fc-l / 1 ' 



c 



if — » oo 



for finite mt 



(3) 



Equation ^ shows that the difference at 
high energy between the exact and approximate 
behaviour is larger at higher orders, so one 
might expect the relative accuracy of the in- 
finite mt approximation to become accordingly 
worse. In Ref. [1^ we have recently computed 
the leading high energy logarithms Eq. ([3]) at 
finite mt, using the techniques of high-energy 
(or kx) factorization [8] (the corresponding co- 
efficients in the mt — *■ oo limit had been previ- 
ously computed in Ref. 0). We can use this re- 
sult to construct an improvement of the NNLO 
result [Sj, by replacing its spurious double loga- 
rithmic growth with the correct high energy be- 
haviour, Eq. Q. 

This construction requires a suitable match- 
ing procedure, and it gives us an approximation 
to the exact NNLO result. We shall first per- 
form this improvement on the NLO contribution, 



Figure 2. Same as Fig. [1] but with 
mH = 280 GeV (here tq 0.315 and lu = 1/20). 



where the exact result is known: it turns out to 
give an approximation to the exact NLO partonic 
cross section which is everywhere accurate to bet- 
ter than 1 % for any value of the Higgs mass, thus 
leading to an approximation to the total cross 
section which is accurate to the level of 0.05%. 
We shall than construct a similar improvement 
of the NNLO term. This improvement changes 
the total cross section computed up to NNLO by 
an amount which varies between 0.2% for light 
Higgs {mn ~ 130 GeV) to 1% for heavy Higgs 
{mH ~ 280 GeV). This is thus the size of the 
error which is made if the cross section is com- 
puted at NNLO using the approximate mt — > oo 
result. By varying the matching prescription, we 
estimate that the ambiguity on this result is at 
the level of the per mille. 

The perturbative expansion of the coefficient 
function in the gluon-gluon channel is 



C(a,;r,2/t)=<5(l-T) + ^CW(T,2/t) 



+ (^) C(^'(r,,t) + 



O (al) 



(4) 



The leading high energy behaviour of the NLO 
and NNLO contributions is given by 

C^'\T,yt)=c['\yt)CA + 0{T) (5) 
C(2)(T,2/t) = -Ci^\yt)Cl\nT + c['\yt) + 0{T). 
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Figure 3. The NNLO coefficient function for 
tuh ~ 130 GeV. The black curve corrdesponds 
to mt — s- oo and the blue one to the approxima- 
tion Eq. dS]) with To = 0.011 and w = 1/100. The 
exact case is not known at NNLO. 



f 2) 

The leading NNLO coefficient C^2 (Vt) 

was com- 
puted in Ref. [TUl, while the subleading NNLO 

(2) 

coefficient (j/t) is unknown. 

The approximate pointlike determination of 
the coefficient function can be improved by sub- 
tracting its spurious small r growth and replacing 
it with the exact behaviour: 



C^PP- (r, yO =C(r,^)+ T(t, to) X 



C(t, yt) - lim Cq{t, qo) 

r— »0 



(6) 



where Co(t, oo) is the sum of contributions to the 
infinite result C{t, oo) which do not vanish 
when T ^ 0, as given at NLO and NNLO by 
the terms listed in Eq. ([5]). Also, T(t, tq) is a 
matching function, which is introduced in order 
to tune the point tq where the small r behaviour 
sets in. We choose 



1 + tanh 



To - T 



(7) 



which in the limit of vanishing width uj becomes 
the step function: lim^^Q T{t, tq) = 8(t — tq). 

At NLO, because the exact asymptotic be- 
haviour Eq. ([5]) is the constant c[^\yt), the 



Figure 4. Same as Fig. [3l but with 
ruH 280 GcV (here to = 0.317 and uj = 1/20). 



matching point tq is naturally determined as the 
value of T where the pointlike approximation 
equals this constant. It is clear from Fig. [1] that 
this choice leads to an excellent approximation to 
the exact result: in fact, it is accurate to better 
than 1% for all t. 

At NNLO the exact asymptotic behaviour 
Eq. ([5]) is a linear rise in Iut. Hence, reasoning 
as at NLO, we are led to choose tq as the point 
where the log derivative of the mt —>■ oo curve 

(2)/ 



matches the asymptotic value C2 {y-, 



din 



-C(2)(t,, 



-9C(2)(y,). 



(8) 



This does not fix completely the approximate re- 
sult at NNLO however, because the subleading 

(2) 

constant C\ {yt) is unknown. We fix it by requir- 
ing that the approximate curve Eq. ([6]) be contin- 
uous at T = To even when the matching function 
r(r, To) = 0(t - Tq). The NNLO approximation 
determined thus is compared to the exact result in 
Figs. [31 m for heavy and fight Higgs, respectively. 
A more conservative matching might consist in- 
stead of taking for To the value found at NLO, 
and then determining again the subleading con- 
stant by continuity. The result found in this way 
is actually very close to the previous one, and in 
fact indistinguishable from it in the case of Fig. 21 
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ruH = 130 GeV 


pointlike 
exact 
appr. 


1.800 
1.797 
1.796 


2.140 
n.a. 
2.136 


rriH = 280 GeV 


pointlike 
exact 
appr. 


1.976 
1.958 
1.959 


2.420 
n.a. 
2.394 



Table 1 

The NLO and NNLO K factors Eq. dH), com- 
puted with center-of-mass energy s = 14 TeV. 



the improved approximate NNLO reduces the er- 
ror due to finite mass terms at NNLO to the per 
mille level. This may be relevant in view of re- 
cent progress on the computation of electroweak 
corrections to this process [12] . 

Finally, we observe that less inclusive quanti- 
ties which depend on the r shape of the par- 
tonic cross section can be rather more sensitive 
to finite-mass effects, in particular if they probe 
the small t tail of a coefficient function. An inter- 
esting case in point, which deserves further inves- 
tigation, is the Higgs rapidity distribution [13j . 



Let us now turn to the inclusive cross section 
We define a K factor 

_ crgg{Th,yt,mjj) 



K{Th;yt) 



gg 



(9) 



where cr°g is the LO cross-section Eq. ([T]), com- 
puted with LO parton distributions and LO cou- 
pling constant. The value of the NLO and NNLO 
K factors, determined using the MRST2002 [11] 
gluon distribution in Eq. ([1]) are given in Table (TJ 
at the LHC center-of-mass energy s = 14 TeV. 
In the table the pointlike, exact and approximate 
cases are shown. 

At NLO the discrepancy between the infinite 
top mass approximation and the exact result is 
tiny, less than 1% even for a fairly heavy Higgs. If 
the improved (approximate) NLO result is used, 
this discrepancy is reduced by a factor three. 
At NNLO the inclusion of the correct small r 
dependence of the partonic coefficient function 
changes the K factor by an amount which varies 
between 0.3% for mn = 130 GeV and 1% for 
rriH = 280 GeV. If we modify the matching pre- 
scription by using the NLO value of tq also at 
NNLO the approximate NNLO results of table 1 
change by 0.1 %. We can take this as the error 
which is made by use of the infinite mt NNLO for- 
mula (there is also a dependence on mt in the con- 
tribution to C^^^ which is proportional to 5(1 — r), 
but at NNLO this contribution is relatively small, 
unlike at NLO). Dominant uncertainties on the 
total Higgs cros-section are typically at the per- 
cent level [7] . Varying to in the matching function 
Eq. O between 1/20 and 1/100 the NNLO resuhs 
change by about 0.1 %. We conclude that use of 
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